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Abstract 

The  success  in  using  VHP  and  UHF  frequency  systems  tor  sounding  polar  ice 
sheets  has  been  tempered  by  an  uncertainly  in  the  in-siiu  dielectric  constant  e; 
which  controls  the  effective  velocity  V9  of  an  electromagnetic  wove  propagefling 
in  an  air-ice  mixture.  An  empirical  equation  tor  determining  e'  vs.  density 
(specific  gravity,  p )  were  proposed  in  1 968  by  Robin  etal.  where  e'=  (1  +  0.851 
p)2.  However,  this  expression  has  met  with  uncertainty  because  wide-angle 
radar  retraction  sounding  techniques  have  produced  values  of  e'  that  are  lower 
than  Robin's  equation  predicts.  This  report  discusses  radar  soundings  mode  on 
the  McMurdo  Ice  Shelf,  Antarctica,  and  compares  the  resulting  e'  values  with 
Robin's  equation,  laboratory  measurements  on  Am  and  ice  and  other  expressions 
given  in  the  literature  tor  determining  e'  vs.  the  specific  gravity  of  dry  fim  and  ice. 
Our  findings  indicate  that  the  form  of  Robin's  equation  is  valid.  However,  our 
analysis  also  indicates  the  expression  could  be  slightly  improved  to  read  e'  = 
(1  +  0.845  p)2.  Reasons  are  suggested  as  to  why  previous  wide-angle  radar 
sounding  studies  did  not  reproduce  Robin's  findings. 


Cover:  Radar  transceiver  antenna  being  towed  from  the  front  of  a  boom 
mounted  to  a  tracked  vehicle.  Behind  the  antenna  box  is  a  disk  and 
support  leg  used  to  elevate  the  boom  above  the  snow  surface.  This 
antenna  arrangement  was  used  In  1974  for  profiling  the  brine  layer  in 
the  McMurdo  Ice  Shelf  and  for  crevasse  detection  during  a  transverse  in 
the  Pensacola  Mountains,  Antaractla,  in  search  of  blue  ice  runway  sites. 
The  radar  console  and  recording  equipment  were  located  in  the  vehicle. 
(Photo  by  A.  Kovacs.) 
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A  Reassessment  of  the  In-Situ  Dielectric  Constant  of  Polar  Fim 

AUSTIN  KOVACS,  ANTHON  Y  J.  GOW  AND  REXFORD  M.  MOREY 


INTRODUCTION 

Since  1948,  when  Steenson  (1951)  first  used  a 
100-MHz  pulse  radio  sounding  system  to  estimate 
ice  thickness  on  the  Seward  Glacier  in  Alaska,  "ra- 
diogladology"  has  been  extensively  used  for  mea¬ 
suring  the  depth  of  polar  ice  sheets.  Success  in  using 
sounding  systems  operating  at  VHF  and  UHF  fre¬ 
quencies  stems  in  part  from  the  generally  very  fa¬ 
vorable  electromagnetic  properties  of  polar  fim  and 
ice,  which  allow  for  deep  penetration  at  these  fre¬ 
quencies. 

The  relative  complex  permittivity  e*  of  an  air-ice 
mixture  can  be  described  by  the  Debye  expression 


where  e'  =  the  relative  real  part  of  e*  (or  the  dielec¬ 
tric  constant) 

e"  =  the  loss  factor  or  imaginary  part  of  e* 
Ejs  =  the  relative  static  dielectric  constant 
x  =  the  relaxation  time 
to  =  the  angular  frequency. 


and  |i'  =  real  part  of  the  magnetic  permeability 
Eg  =  real  effective  (bulk)  dielectric  constant 
oe  =  effective  (bulk)  conductivity. 

Since  the  conductivity  of  dry  polar  fim  and  ice  is 
extremely  low,  the  effective  phase  velocity  of  high- 
frequency  electromagnetic  waves  can  be  deter¬ 
mined  from 

Ve-c/^e  (2) 

where  c  is  the  velocity  of  an  electromagnetic  wave  in 
a  vacuum,  0.3  m/ns. 

If  the  two-way  vertical  travel  time  t  of  a  transmit¬ 
ted  wavelet  traveling  from  the  surface  to  the  ice  bot¬ 
tom  (or  from  some  internal  horizon)  and  back  is 
measured,  then  this  time  may  be  used  to  estimate 
the  distance  D  to  the  reflective  boundary  from 

D=_fe=  =  iVe.  (3) 

zVe;  2 

To  accurately  make  this  distance  determination,  e' 
for  the  site  must  be  known.  Eg  can  be  determined 
when  f  and  D  are  known  as  follows: 


In  this  report  the  variation  in  the  dielectric  constant 
vs.  fim  density  or  specific  gravity  is  discussed,  as 
this  value  governs  the  propagation  speed  of  an 
electromagnetic  wave  transmitted  from  VHF  and 
UHF  sounding  systems  into  polar  fim  and  ice.  The 
loss  tangent,  tan  8,  is  represented  by  e"  /e'  and  for 
polar  fim  and  ice  is  typically  »  0.1 

An  electromagnetic  wave  propagates  through 
fim  and  ice  with  an  effective  phase  velocity  Ve  of 

Ve  =  co/p 

where  p  is  the  phase  constant  determined  by 


p  =  © 


H'C1/2 


1+- 


<*e 


\l/2 


or2  (e;r 


+ 1 


1/2 


Since  D  is  seldom  known,  this  equation  is  of  limited 
use.  Where  fim  and  ice  density  profile  data  are 
available  for  the  site,  it  is  possible  to  estimate  e'  for  a 
given  depth  increment  based  on  the  mean  interval 
density  and  a  laboratory-determined  Eg  value  for 
this  mean  density. 

To  circumvent  these  requirements,  borehole  in¬ 
terferometry  and  wide  angle  radar  reflection  tech¬ 
niques  have  been  used  to  determine  D  or  Eg.  These 
methods  have  been  widely  discussed  by  Bentley 
(1979),  Jezek  and  Roeloffe  (1983),  Bogorodsky  et  al. 
(1985)  and  others.  However,  Jezek  and  Roeloffe 
(1983)  state  that  “collectively,  these  experiments  are 
inconclusive."  It  has  been  found  that  Ve  determined 
from  wide-angle  reflection  sounding  measurements 


tends  to  be  too  high,  and  the  resulting  low  value  for 
e'  is  not  in  agreement  with  laboratory  determina¬ 
tions  of  e'  made  on  fim  and  ice  of  the  same  density. 
System  timing  error  and  curved  ray  paths  are  two  of 
the  possible  reasons  given  for  the  observed  discrep¬ 
ancies  (Jezek  and  Roeloffe  1983). 

The  lack  of  agreement  between  laboratory-deter¬ 
mined  e'  values  and  the  various  radar  field  results 
has  brought  into  question  the  validity  of  Robin's 
(Robin  et  al.  1969,  Robin  1975)  borehole  radar  inter¬ 
ferometry  measurements,  from  which  he  proposed 
an  empirical  equation  relating  the  index  of  refrac¬ 
tion  n  to  the  fxm-ice  density,  which  is  henceforth  ref¬ 
erenced  as  the  specific  gravity  p  for  the  purpose  of 
making  his  and  subsequent  equations  dimension- 
ally  correct: 


n  =-£-  =  1  +  0.851  p. 

V'e 

Since  n2  =  t'e  in  fire  radio-radar  frequency  range  of 
interest,  then 

e;  =  (l +0.851  p)2.  (5) 

The  form  of  this  exponential  expression  is  hereafter 
referred  to  as  a  refraction  typ  xjuation. 

During  their  study  of  the  -me  infiltration  zone 
in  the  McMurdo  ke  Shelf,  Antarctica  (Fig  1),  Kovacs 
and  Gow  (1975)  and  Kovacs  et  al.  (1981, 1982)  used 
two  impulse  radar  sounding  methods  to  profile  fire 
depth  variations,  lateral  continuity  and  "inland" 
boundary  of  seawater  infiltration  into  fire  ice  shelf 


Figure  1.  Map  of  the  McMurdo  Ice  Shelf  area.  The  ice/fim  transition  controls  the  western  limit  of  brine 
infiltration.  The  eastern  limit  of  brine  infiltration  is  controlled  by  ice  porosity  and  westward  ice  shelf  move¬ 
ment.  Very  close  to  the  western  side  of  White  Island  snow  accumulation  can  exceed  2  m/yr  and  the  ice  shelf  is 
composed  of  recently  infiltrated  snow  ice.  In  the  crevassed  area  at  the  north  end  of  White  Island  is  a  small 
population  of  land-locked  seals  (from  Kovacs  et  al.  1981). 
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Surface 


Figure  3.  Radar  records  ofMcMurdo  Ice  Shelf  along  sounding  line.  Top  graphic  record  shows  300-MHz  radar  profile  of  the  brine  layer  along  the  sounding  line  shown  in 
Figure  1.  Steps  in  the  brine  layer  are  the  result  of  ice  shelf  bredcout  events  that  led  to  a  new  zone  of  seawater-infiltrated  fim  abooe  the  older  brine  layer.  Note  the  internal  fim 
layering,  especially  on  the  right  side  of  the  graphic  record.  The  bottom  record  shows  the  terminus  of  the  brine  zone  and  the  sloping  bottom  of  the  McMurdo  Ice  Shelf.  (From 
Kovacs  et  al.  1981.) 


(Fig.  2).  Ok  method  previously  discussed  (Kovacs 
et  aL  1982)  used  a  fixed  separation,  dual-antenna 
configuration  in  which  one  antenna  was  a  trans¬ 
ceiver  (colocated  transmitter-receiver  antenna)  and 
the  second  antenna  was  a  receiver  only.  The  dual¬ 
antenna  methodology  used  is  discussed  in  Kovacs 
and  Morey  (1979).  YTith  this  radar  sounding  meth¬ 
od,  Ve,  D  and  were  determined  for  tire  fim  be¬ 
tween  the  snow  surface  and  the  brine  infiltration 
layer  The  depth  of  the  brine  layer  was  determined 
by  direct  borehole  measurements  one  year  after 
the  radar  measurements  were  made.  After  fire  new 
snow  accumulation  was  taken  into  consideration, 
the  agreement  between  fire  radar-  and  borehole- 
determined  depths  was  typically  within  2%.  These 
results  indicate  that,  for  the  shelf  fim  layer,  where 
the  most  significant  depth  variation  in  density  oc¬ 
curs,  the  dual-antenna  radar  sounding  method  pro¬ 
vided  a  satisfactory  means  for  determining  t'e  and 
Ve  as  well  as  the  depths  to  internal  layers,  which 
were  dearly  revealed  in  the  radar  records  (Fig.  3 
and  4). 

While  Kovacs  et  aL  (1982)  showed  good  agree¬ 
ment  between  their  radar-determined  depth  to  die 
brine  layer  and  the  drill-hole-measured  depth  and 
gave  estimates  of  e'  and  Ve,  they  did  not  provide  in¬ 
formation  for  an  assessment  of  the  Robin  equation. 
Nor  did  they  attempt  "to  resolve  the  question  of  the 
possible  discrepancy  between  laboratory  and  field 
dielectric  constant  measurements"  for  fim  and  ice  as 
did  Jezek  and  Roeloffs  (1983).  In  this  report,  data  ob¬ 
tained  from  the  vertical  two-way  impulse  radar 
sounding  measurements  on  the  McMurdo  Ice  Shelf 
are  presented.  The  findings  are  discussed  in  relation 
to  the  Robin  equation  and  other  field,  laboratory, 
and  theoretical  determinations  of  the  dielectric  con¬ 
stant  of  an  air-ice  mixture. 


FIELD  MEASUREMENTS 

A  Geophysical  Survey  Systems  impulse  radar 
system  was  used  in  the  field  study.  The  transceiver 
antenna,  when  snow-coupled,  transmitted  a  broad¬ 
band  wavelet  with  a  carter  frequency  of  about  100 
MHz.  The  vertical  two-way  travel  time  of  the  wave¬ 
let  from  the  snow  surface  to  various  subsurface  an¬ 
nual  accumulation  layers  and  tire  top  of  the  brine 
infiltration  zone  in  the  McMurdo  Ice  Shelf  was  re¬ 
corded  in  real  time  on  a  graphic  recorder  (Fig.  4)  and 
on  a  magnetic  tape  recorder  for  later  analysis  on  a 
frequency  analyzer.  The  two-way  flight  time  was 
measured  from  where  tire  transmit  wavelet  first 
crossed  the  zero  voltage  line  to  where  tire  relected 
wavelet  of  interest  first  crossed  tire  zero  line  as  de¬ 
picted  in  Figure  5. 

Two-way  travel  time  measurements  were  made 
at  six  sites  in  1978  where  cores  were  obtained  for  de¬ 
termining  tire  depth  to  the  brine  layer  and  the 
depth-density  profile.  The  measured  distance  from 
the  edge  of  tire  McMurdo  Ice  Shelf,  brine  layer 
depth,  average  fim  density  and  the  two-way  travel 
time  are  listed  in  Table  1  as  are  the  calculated  Eg  and 
Ve  values  for  each  radar  sounding  site.  Also  given  is 
the  location  and  the  two-way  travel  time  to  tire  bot¬ 
tom  of  tire  ice  shelf  at  station  307,  which  was  located 
immediately  beyond  tire  1978  brine  infiltration  limit 
(Figs.  2  and  3).  The  depth  of  the  ice  shelf  at  station 
24,  located  just  40  m  before  station  307,  was  esti¬ 
mated  by  Kovacs  et  aL  (1982).  Using  the  1977  dual- 
antenna  radar  sounding  data,  they  first  determined 
the  depth  to  tire  brine  layer  at  station  24.  The  two- 
way  flight  time  from  tire  brine  layer  to  the  bottom  of 
the  ke  shelf  was  then  used,  along  with  an  assumed 
value  of  Eg  =  2.95  for  this  depth  interval  The  thick¬ 
ness  of  tire  ice  shelf  was  then  estimated  (using  eq  2 


Table  1. 1978  McMurdo  Ice  Shelf  station  data  along  traverse  from  shelf  edge  to  ice  shelf  move¬ 
ment  marker  307. 


Station 

no. 

Distance 

inland 

(m) 

Brine 

depth* 

(m) 

Average 
specific 
gravity  p 

Wavelet 
travel  time  t 
(ns) 

Effective  bulk 
dielectric 
constant  r£ 

Effective 
velocity  Vf 
(m/ns) 

A 

5 

23 

0.395 

22 

1.74 

0.227 

B 

770 

835 

0509 

83 

1.98 

0.213 

10b 

2900 

13  3 

0540 

130 

2.08 

0.208 

C 

3590 

9.45 

0565 

190 

2.15 

0.204 

D 

7370 

33.8 

0.624 

344 

2.33 

0.196 

E 

9560 

50.4 

0.686 

529 

2.48 

0.190 

307 

9700 

— 

— 

1263 

— 

— 

’Drillhole  measurements 


5 


Figure  4.  Graphic  record  of 300-MHz  radar  sounding  data  showing  internal  layers  in  the  snow  field  below  Castle  Rock 
(Fig.  J ).  The  dotted  line  marks  the  site  where  a  20-rn  firn  core  was  obtained.  Core  analysis  revealed  the  internal  layers 
were  summer  stratigraphy  horizons. 
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Figure  5.  Example  of  recemsd  radio  echo 
wavelet.  In  the  graphic  record  many  horizons 
or  internal  layers  were  observed  in  the  ice  shdf 
fim.  These  layets  could  be  tracked  for  many 
kilometers.  Their  voltage  amplitude  signature 
is  shown  above  (from  Kooacs  et  al.  1982). 


Table  2.  McMurdo  Ice  Shelf  depth  increment  data. 


Station 

increment 

Depth 

increment 

(m) 

Average 
specific 
gravity  p 

Wavelet 
travel  timet, 
(ns) 

Effective  bulk 
dielectric 
constant 

Effective 
velocity  V, 
(m/ns) 

A-B 

635 

0304 

61 

2.08 

0206 

B-C 

10.6 

0.633 

107 

229 

0.199 

C-D 

1435 

0.726 

154 

239 

0.186 

D-E 

16.6 

0.808 

185 

2.79 

0.180 

and  3)  to  be  113  m.  This  thickness 
estimate  will  be  compared  with 
the  estimated  thickness  for  sta¬ 
tion  307  as  determined  from  the 
1978  radar  two-way  travel  time 
data  obtained  at  this  site. 

The  brine  layer  depths  and 
two-way  travel  time  differences 
between  stations  A-B,  B-C,  C-D 
and  D-E  are  listed  in  Table  2 
along  with  foe  average  fim  densities  for  the  depth 
increments.  Also  listed  are  foe  calculated  e'  and  Ve 
values  for  the  depth  intervals.  A  plot  of  the  e'  vs.  p 
values  given  in  Tables  1  and  2  is  shown  in  Figure  2. 

Numerous  investigators  have  determined  that,  at 
— 10°C  and  a  p  of  0.917,  the  dielectric  constant  of 
pure  polycrystalline  ice,  at  foe  frequencies  of  inter¬ 
est,  is  well  represented  by  a  value  of  3.15  (e.g., 
Cummings  1952,  Evans  1965,  Ulaby  et  al.  1986,  Koh 
1992).  Koh  (1992)  also  reaffirmed  that  e'  is  high-fre¬ 
quency-independent.  Fujita  et  al.  (1992)  reconfirmed 
that  e'  for  polycrystalline  ice  is  also  insensitive  to 
temperature,  varying  at  about 0.0009 /C°.  Therefore, 
in  foe  subsequent  discussion  and  equations,  an  e' 
value  of  3.15  will  be  used  for  pure  bubble-free  ice  at 
a  specific  gravity  of  0.917.  The  intercept  of  these  val¬ 
ues  is  shown  by  a  cross  in  Figure  6  and  in  subse¬ 
quent  figures. 

It  has  been  found  that  when  glacial  ice  is  under 
stress,  as  in  a  flow  regime,  an  anisotropy  develops  in 
which  the  c-axes  of  the  ice  crystals  become  aligned. 
Johari  and  Chasette  (1975)  found  that  in  the  infrared 
region  a  dielectric  constant  change  on  the  order  of 
0.004  could  arise  if  the  electric  field  was  parallel 


rather  titan  perpendicular  to  the  crystal  c-axis.  Simi¬ 
larly,  Fujita  et  al  (1992)  have  shown  experimentally 
that  the  change  could  be  as  much  as  0.04  at  9.7  GHz. 
Thin-section  studies  of  the  crystalline  structure  of 
the  cores  obtained  from  foe  McMurdo  Ice  Shelf  re¬ 
vealed  that  the  ice  grains  had  a  random  c-axis  fabric 
that  existed  from  the  surface  to  foe  bottom  of  our  50- 
m-deep  drill  hole  near  station  307  (Fig.  2).  This  ran¬ 
dom  fabric  orientation  therefore  eliminates  ice  crys¬ 
tal  anisotropy  effects  as  a  contributing  factor  in  our 
radar  measurements. 

Three  curves  are  also  shown  in  Figure  6.  One 
curve  represents  Robin's  equation  (eq  5).  Another  is 
a  refraction  type  curve  statistically  fitted  through  the 
data.  The  equation  for  this  curve, 

e'r =(0.992  +  0.848  p)2  (6) 

fits  foe  data  with  an  r2  =  0.989  and  a  standard  error 
of  0.035.  The  third  curve  shown  was  determined  as 
follows.  The  four-parameter  volume  fraction  type 
equation 
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p.  Specific  Gravity 


Figure  6.  Field-determined  insitu  effective 
dielectric  constant  e’e  vs.  the  bulk  specific 
gravity  ofMcMurdo  Ice  Shelf fim.  The  curves 
passing  through  the  field  data  •  were  derived 
from  the  indicated  equations  which  are  described 
in  this  report.  The  +  represents  the  dielectric  con¬ 
stant  (3.15)  of  pure  ice(p  =  0.917). 


where  e„  is  the  dielectric  constant  of  air  (=  l)ande'is 
the  dielectric  constant  of  pure  ice,  va  and  die  vol¬ 
ume  fractions  of  air  and  ice,  respectively,  were  used 
to  calculate  the  dielectric  constant  of  fim  and  ice 
over  die  specific  gravity  range  from  0  to  0.917  in  in¬ 
crements  of  05.  A  refraction  type  expression  was 
that  statistically  fitted  through  this  data.  From  this 
procedure  it  was  found  that  eq  7  is  well  represented, 
at  a  cross-correlation  of  r2  =  1.000,  by  the  expression 

c'  =  (l  +  0.845p)2.  (8) 

It  is  apparent  from  Figure  6  that  both  Robin's  em¬ 
pirical  equation  (eq  5)  and  the  simplified  volume 
fraction  expression  (eq  8)  fit  the  McMurdo  Ice  Shelf 
dielectric  constant  vs.  specific  gravity  field  data  very 
well. 

Numerous  other  empirical  and  theoretical  ex¬ 
pressions  have  been  used  to  describe  the  dielectric 
constant  of  an  air-ice  mixture.  A  number  of  these  ex¬ 
pressions  will  be  compared  with  the  McMurdo  Ice 
Shelf  data  and,  where  statistically  permissible,  die 
expressions  will  be  simplified  to  a  refraction  equa¬ 
tion,  following  the  procedure  used  to  obtain  eq  8. 

Looyenga  (1965)  proposed  a  mixing  formula 
based  on  spherical  inclusions 

er  =  fcr'i 73  -  eJ73)  Vi  +  el73]3  (9) 

and  because  Ea  =  1  then 

er  =  [(eri 73  -  l)  vi  +  l]3- 


e'r  =  (1  +  0.508 p)3  (10) 

Equation  9  can  also  be  well  represented  by  a  refrac¬ 
tion  type  equation,  with  an  r2  =  0.999,  as 

e;  =  (0.988  +  0.845 p)2  (11) 

which  is  in  very  near  agreement  with  eq  8. 

Wiener  (1910)  presented  a  two-phase  mixing  for¬ 
mula  containing  an  empirical  structure  factor  called 
die  Formzahl  p.  This  function  varies  with  the  shape 
and  orientation  of  the  inclusions.  In  fim  and  ice  the 
Formzahl  may  vary  if  the  medium  is  nonisotropic. 
The  form  of  the  equation  for  an  air-ice  mixture  is 

Er  ~  1  =  Vj  (  E»—  | 

e'+p  \e-+p/ 

+  (l-Vi)pJl!).  (12) 

Since  Eg  is  one,  the  second  term  on  the  right  is  zero.  It 
was  found  through  statistical  regression  that  eq  12 
fits  the  McMurdo  Ice  Shelf  fim  dielectric  constant  vs. 
specific  gravity  data  best  when  p  =  5.5.  This  value  is 
in  good  agreement  with  die  values  of  5  given  in 
Bogorodsky  et  al.  (1985)  and  75  given  in  Ambach 
and  Denoth  (1972)  for  vertically  propagating  elec¬ 
tromagnetic  waves  in  fim  and  ice.  But  it  is  not  in 
agreement  with  Shabtaie  and  Bentiey  (1982),  who 
found  that  p  =  0  provided  the  best  fit  for  their  com¬ 
mon  reflection  point  measurements  at  Dome  C,  Ant¬ 
arctica. 

The  results  from  this  measurement  technique 
have  caused  the  most  concern  regarding  die  validity 


This  equation  was  simplified  by  Stiles  and  Ulaby 
(1981)  to  read 
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Figure  7.  Field-determined  in-situ  effective  di¬ 
electric  constant  e'  os.  the  bulk  specific  gravity 
of  McMurdo  Ice  Shelf  fim.  The  curves  passing 
through  the  field  data  were  derived  from  the  indi¬ 
cated  equations  which  are  described  in  this  report. 


of  the  Robin  equation.  It  is  generally  expected  that 
rpolar  fim  and  ice  p.  should  be  above  2  (e.g.,  Evans 
1965).  A  Formzahl  in  the  range  of-  5-7 suggests  that 
the  horizontal  layering  in  polar  fim  has  a  weak  ef¬ 
fect  on  e' .  Equation  12  can  be  simplified  to 

e'  =  (1  +  0.840p)2  (13) 

which  cross-correlates  extremely  well  with  the 
Wiener  equation,  when  it  =  5.5,  at  t2- 1.000. 

The  curves  for  eq  10, 11  and  13  and  the  McMurdo 
Ice  Shelf  e'e  vs.  p  data  are  shown  in  Figure  7.  Here 
again  the  curves  for  the  simplified  equations  are 
seen  to  fit  the  field  data  very  well. 

A  mixing  formula  by  deLoor  (1956),  when  ap¬ 
plied  to  an  air-ice  medium  containing  spherical  in¬ 
clusions,  is 


=  va 


ea~l 

£a  +  2e' 


+  Vi 


1 

le'j  +  2e' J 


Since  £3  =  1,  the  first  term  on  the  right  side  is  zero, 
and  the  equation  for  an  air-ice  mixture  becomes 


(16) 


which  in  turn  is  well  represented  at  a  cross-correla¬ 
tion  of  r2  =  0.999  by 

e'  =  (0.988  +  0.850p)2.  (17) 

For  disk-shaped  particles  the  Polder  and  van 
Santen  (1946)  equation  is 


[3(1-4)  (e'n-D-e'ri 


(14) 


c>  -  3gfi  +  2e'j  (e',  - 1)  (1-40 

e'ri-(e'ri-l)(l-<») 


+ 


2±V[3(l-<tO  (e;i-l)-e;i+2]2 


which  can  be  simplified  to 


where  $  is  the  porosity.  This  equation,  when  simpli¬ 
fied  to  a  refraction  type  expression,  becomes 

e;  =  (0.983  +  0.858 p)2  (15) 

with  a  cross-correlation  of  r2  =  0.999. 

The  Bottcher  (1933)  and  Polder  and  van  Santen 
(1946)  mixing  formula  for  a  medium  composed  of 
spherical  inclusions  is 


e;  +  (1.006  +  0.839 p)2  (18) 

with  a  cross-correlation  of  r2  =  1.000. 

The  Tinga  et  aL  (1973)  mixing  formula  for  an  air- 
ice  medium  containing  spherical  inclusions  is 


e'r  =  1  + 


avife-i) 


(2  +  3e'i)-v,(e'i-l) 


(19) 
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Figure  8.  Field-determined  in-situ  effective  di¬ 
electric  constant  e'us.  the  bulk  specific  gravity  of 
McMurdo  Ice  Shelf  fim.  The  curves  passing 
through  the  field  data  were  derived  from  the  indi¬ 
cated  equations  which  are  described  in  this  report. 


which  when  simplified  becomes 


which  can  be  simplified  to 


£'r  =  (1  +  0.910  p)  /(.l  -  0.455  p)  (20) 


e;  =  0.367  +  3.035  p 


(24) 


with  a  cross-correlation  of  r2  =  1.000.  Equation  19 
does  not  cross-correlate  very  well  with  a  refraction 
type  expression. 

The  curves  respresenting  eq  15,  17  and  20  are 
shown  in  Figure  8,  along  with  the  e'e  vs.  p  data  from 
the  McMurdo  Ice  Shelf.  Here  we  see  that  the  curves 
representing  eq  15  and  17  pass  nicely  through  the 
field  data  but  the  curve  representing  eq  20  does  not. 
From  Figure  8  it  appears  that  the  Tinga  et  al.  equa¬ 
tion  should  not  be  used.  Sihvola  et  al.  (1985)  also 
found  that  the  Tinga  et  al.  expressbn  (eq  19)  is  not 
appropriate  for  determining  e'  for  an  air-ice  mix¬ 
ture. 

There  are  other  expressions  that  also  do  not  fit  the 
McMurdo  Ice  Shelf  field  data  well.  Three  of  these 
equations  are  given  below.  The  first  is  a  volume  av¬ 
erage  expression  (Lang  1983): 

Er  =  $ea  +  (l  -^)£ri  (21) 

which  can  be  further  simplified  to 

e'r=l+  2.345  p  (22) 

with  a  cross-correlation  of  r2  =  1.000. 

Ackley  and  Keliher  (1979)  proposed  the  use  of  a 
mixture  formula  by  van  Beek  to  determine  e'  for  po¬ 
lar  ice  containing  spherical  inclusions: 


£r  =  «ri 


1  + 


^(Ea-E'ri) 
£ri  +  £a 


(23) 


with  a  cross-correlation  of  r2  =  1.000. 

The  linear  model  proposed  by  Hufford  (1991)  for 
an  air-ice  mixture  is 


e;  =  l  +  3Vi  p— M  (25) 

\eri  +  2/ 

which  can  be  reduced  to 

e'  =  1  +  1.366p  (26) 

with  a  cross-correlation  of  r2  =  1.000. 

Equations  22, 24  and  26  are  graphically  shown  in 
Figure  9,  along  with  eq  6  and  the  McMurdo  Ice  Shelf 
Eg  vs.  p  field  data.  Clearly,  eq  22  does  not  track 
through  the  field  results  and  should  not  be  used. 
Equations  24  and  26  appear  to  have  a  limited  range 
of  use;  however,  they  are  not  considered  acceptable 
in  light  of  how  well,  for  example,  eq  6  fits  the  data. 

Other  expressions  for  determining  the  dielectric 
constant  of  an  air-ice  mixture  vs.  specific  gravity  are 
listed  in  Appendix  A. 

To  further  show  that  a  simple  refraction  type 
equation  can  well  represent  the  dielectric  constant  of 
a  mixture  of  air  and  ice,  the  10  McMurdo  Ice  Shelf 
dielectric  constant  vs.  specific  gravity  data  from 
Tables  1  and  2  were  combined  with  15  laboratory 
determinations  of  Cummings  (1952)  at  9.4  GHz,  17 
field  measurements  made  at  20  MHz  on  "old"  al¬ 
pine  snow  by  Denoth  (1978)  and  47  laboratory  mea¬ 
surements  by  Nyfors  (1982)  made  at  UHF  and  SHF 
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Figure  9.  Field-determined  in-situ  effective  di¬ 
electric  constant  e'  vs.  the  bulk  specific  gravity 
ofMcMurdo  Ice  Shelf  Jim.  The  curves  passing 
through  the  field  data  were  derived  from  the  indi¬ 
cated  equations  which  are  described  in  this  report. 


Figure  10.  Field  and  laboratory-determined 
values  of  the  dielectric  constant  of fim  vs.  spe¬ 
cific  gravity.  The  regression  curve  passing 
through  the  data  is  represented  by  equation  27  as 
discussed  in  this  report. 


frequencies  on  a  variety  of  snow  types  (e.g.,  fine- 
and  coarse-grained  snow).  A  plot  of  these  data 
along  with  a  weighted  point  ofe'j  =  3.15  at  p  =  0.917 
is  given  in  Figure  10.  The  resulting  regression  curve 
passing  through  the  data  is  derived  from 

e'  =  (1.000  +  0.845 p)2  (27) 

This  curve  fits  the  data  with  a  correlation  of  r2  = 
0.999  and  a  standard  error  of  ±  0.031.  It  should  be 
noted  that  eq  27  is  also  the  simplified  form  of  the 
volume  fraction  expression  (eq  8). 

The  justification  for  adding  a  weighted  dielectric 
constant  of  3.15  at  a  specific  gravity  of  0.917  to  the 
data  in  Figure  7  is  based  on  our  earlier  determina¬ 
tion  that  this  dielectric  constant  is  a  representative 
value  for  pure  ice.  In  addition,  the  weighted  value 
was  added  because  of  the  paucity  of  dielectric  con¬ 


stant  test  data  at  the  higher  specific  gravities.  With¬ 
out  the  weighted  point,  the  regression  equation 
would  read  e';  =  (1.002 + 0.836  p  j2  with  an  r2  of  0.994 
and  a  standard  error  of 0.035.  This  equation  gives  an 
e'  of  3.13  at  a  p  of  0.917,  the  same  value  obtained  by 
using  eq  6. 

It  was  previously  stated  that  Jezek  and  Roeloffs 
(1983)  made  radar  sounding  measurements,  in 
Greenland  and  Antarctica,  for  the  purpose  of  evalu¬ 
ating  the  validity  of  the  Robin  equation.  In  these 
measurements,  a  metal  target  was  lowered  down  a 
borehole  and  tire  travel  time  of  a  transmitted  wave 
from  the  surface  to  the  target  and  back  was  re¬ 
corded.  Using  the  known  depth  to  the  target  and  fire 
density  of  the  fim  and  ice  above  the  target,  they  cal¬ 
culated  the  effective  radar  wave  velocity  in  fire  me¬ 
dium  for  different  target  depths.  The  average  Ve 
data  for  each  depth  at  the  Dye  3  Greenland  site  are 
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Va.  Effective  Vatodly  (n^i») 
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Figure  11.  Effective  radar  velocity  vs.  depth  at  DYE-3 
Greenland.  The  data  points  are  jezek  and  Roeloffs'  average 
values,  the  Robin  curve  is  from  their  analysis  and  the  right 
curve  is  based  oneql7  in  this  report. 


shown  in  Figure  11.  Also  shown  is  the  Vevs.  depth 
curve  they  calculated  using  Robin's  equation,  the 
depth-density  profile  data  for  the  site  and  eq  3.  The 
regression  curve  we  fitted  to  their  velocity  data  was 
selected  because  it  paralleled  the  calculated  curve 
using  Robin's  equation. 

As  indicated  in  Figure  11,  the  regression  curve 
through  Jezek  and  Roeloffs'  data  gives  an  apparent 
2  m/ps  slower  velocity  gradient  than  the  velocity 
vs.  depth  curve  calculated  with  use  of  Robin's  equa¬ 
tion. 

The  third  curve  shown  in  Figure  11  was  calcu¬ 
lated  with  the  use  of  eq  27,  the  site  depth-density 
data  and  eq3.  As  one  would  expect,  this  curve  gives 
a  slightly  higher  velocity  vs.  depth  profile.  Never¬ 
theless,  the  agreement  between  the  three  curves  is 
considered  to  be  very  good. 

We  believe  that  the  close  similarity  between  eq  5 
and  eq  6  and  27  reaffirms  that  the  empirical  assess¬ 
ment  of  Robin  et  al  (1969)  and  Robin  (1975)  for  de¬ 
scribing  the  dielectric  constant  of  polar  fim  and  ice 
vs.  their  specific  gravity  is  valid.  It  may  be  of  interest 


to  note  that  the  form  of  this  refraction  expression 
may  also  be  used  to  describe  the  dielectric  constant 
of  soil  vs.  its  moisture  content,  as  shown  in  Appen¬ 
dix  B. 


DISCUSSION 

The  borehole-measured  depths  to  the  brine  layer 
in  the  McMurdo  Ice  Shelf  were  plotted  vs.  the  radar 
wavelet  two-way  flight  time  measurements  in  Fig¬ 
ure  12.  The  two-way  flight  time  measured  to  the  bot¬ 
tom  of  the  ice  dielf  in  1978  at  station  307  was  1263  ns 
(Table  1).  From  the  regression  equation  given  in  Fig¬ 
ure  12,  this  time  would  represent  a  depth  of  116  m. 
As  previously  stated,  Kovacs  et  al.  (1982)  gave  an 
estimated  thickness  of  113  m  for  the  shelf  at  station 
24  in  1977.  Since  station  307  was  40  m  east  of  station 
307,  one  could  assume  that  the  ice  is  ~  03  m  thicker 
at  station  307  than  at  station  24  because  the  ice  shelf 
increases  in  thickness  to  the  east  of  this  station  (Fig. 
2)  at  about  8  m/km.  In  any  event,  the  depth  differ¬ 
ence  between  the  two  determinations  (113  m  and 
116  m)  is  about  a  21/2%  variance.  Clearly  the  two 
assessments  of  ice  shelf  thickness  are  in  close  agree¬ 
ment,  given  the  snow  accumulation  (0.9  m  between 
surveys),  bottom  melting  and  shelf  movement  dy¬ 
namics,  which  can  cause  short-  as  well  as  long-term 
ice  shelf  thickness  variations  at  a  location  that  is  not 
fixed  in  space.  In  addition,  the  1982  estimate  was 
based  on  an  assumed  bulk  ice  density  for  the  shelf 
below  the  brine  layer  depth,  a  density  that  now  ap¬ 
pears  to  be  a  bit  high. 

Another  interesting  plot  is  that  of  the  dielectric 
constant  values  vs.  depth,  as  shown  in  Figure  13. 
From  the  regression  equation  for  the  curve  in  Figure 
13,  for  the  ~116-m-thick  ice  shelf  at  station  307 can 
be  estimated  to  be  -  2.74.  Jericek  and  Bentley  (1971) 
gave  an  Eg  of  3.01  ±  0.03  for  the  510-m-thick  Ross  Ice 
Shelf  near  Roosevelt  Island,  and  Jezek  et  al.  (1978) 
estimated  an  Eg  value  of  2.99  for  the  480-m-thick 
Ross  Ice  Shelf  at  station  N19.  The  equation  in  Figure 
13  gives  an  Eg  of  3.02  and  3.01  for  these  two  ice  shelf 
thicknesses,  respectively.  The  dose  agreement  be¬ 
tween  the  field  and  calculated  values  are  surprising 
when  one  considers  that  the  depth-density  profile 
of  the  fim  layer  can  be  expected  to  vary  with  ice 
shelf  locations.  Nevertheless,  this  variation  may  not 
be  enough  to  significantly  influence  the  bulk  dielec¬ 
tric  constant  of  the  fim  layer  and  in  turn  that  of  the 
ice  shelf.  If  so,  then  the  results  from  the  McMurdo 
Ice  Shelf,  a  lobe  of  the  Ross  Ice  Shelf,  may  be  appli¬ 
cable  to  other  areas  of  the  Ross  Ice  Shelf.  In  any 
event.  Figure  12  may  be  used  for  estimating  die 
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Figure  12.  Horizon  depth  within  the  Mc- 
Murdo  Ice  Shelf  vs.  radar  wavelet  two-way 
flight  time. 


0,  Depth  (m) 


Figure  13.  Effective  dielectric  constant  of  the 
fim-ice  in  the  McMurdo  Ice  Shelf  vs.  depth  in¬ 
crement. 


depth  of  these  ice  shelves  from  the  two-way  travel 
times  obtained  from  vertical  radar  sounding  mea¬ 
surements. 

However,  a  note  of  caution  is  in  order.  For  ex¬ 
ample,  the  Ross  Ice  Shelf  at  Site  J-9  is  known  to  have 
a  ~6-m-thick  layer  of  saline  ice  accreted  to  the  bot¬ 
tom  (Zotikov  1979,  Zotikov  et  al.  1980,  Morey  and 
Kovacs  1982).  Because  this  ice  attenuates  electro¬ 
magnetic  energy  at  the  frequencies  often  used  in  ra¬ 
dar  sounding  (as  does  the  brine-soaked  fim),  and 
given  the  spreading  losses  which  occur  with  dis¬ 
tance,  a  reflection  from  the  saline  ice/seawater  inter¬ 
face  may  not  be  obtained.  Nevertheless,  a  reflection 
may  be  received  from  the  fresh  ice/saline  ice  inter¬ 
face,  and  therefore  foe  depth  to  this  interface  could 
be  estimated  using  a  two-way  travel  time  measure¬ 
ment  and  Figure  12.  Should  a  reflection  also  be  ob¬ 


tained  from  foe  saline  ice/seawater  interface,  then 
foe  two-way  travel  time  measured  in  this  saline  ice 
layer  may  be  used  to  estimate  foe  layer  thickness  T 
from 

T  =  -0.708  +  0.088 1  (28) 

and  foe  layer  effective  dielectric  constant  from 

e;  =  5.983  exp(-0.477T)  +  3.172  (29) 

Equations  28  and  29  are  from  foe  results  of  impulse 
radar  two-way  time  of  flight  measurements  made  at 
-  100  MHz  on  thick  sea  ice  by  Kovacs  and  Morey 
(1990).  For  saline  ice  growth  over  10  m  thick  on  tire 
bottom  of  an  ice  shelf,  a  value  of  325  ±  0.05  for  e' 
may  be  reasonable  for  estimating  T  using  eq  3 
wherein  T  is  substituted  for  D. 
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It  set  ms  that  many  distributions  in  nature  tend  to 
follow  exponential-like  trends.  The  above  results  for 
the  variation  in  the  dielectric  constant  of  fim  and  ice 
vs.  specific  gravity  as  well  as  the  results  for  the  di¬ 
electric  constant  vs.  the  moisture  content  of  soil 
(App.  B)  appear  to  represent  two  of  them. 

CONCLUSIONS 

The  findings  of  this  report  reaffirm  that  a  simple 
refraction  type  equation  can  be  used  to  describe  ex¬ 
tremely  well  the  dielectric  constant  vs.  specific  grav¬ 
ity  of  fim  and  ice.  It  was  further  verified  that  the  so- 
called  Robin  expression  [e'  =  (1  +  0.851  p)2]  is  very 
suitable  for  determining  the  dielectric  constant  vs. 
specific  gravity  of  dry  polar  fim  and  ice,  but  that  it 
may  be  slightly  improved  by  changing  the  constant 
0.851  in  the  equation  to  0.845,  especially  if  one  ac¬ 
cepts  3.15  for  the  dielectric  constant  of  pure  ice  as 
measured  by  Cummings  (1952)  again  reaffirmed  by 
Koh  (1992).  Why  previous  investigators'  field  data 
did  not  correlate  well  with  the  Robin  equation  is  de¬ 
batable.  Jezek  and  Roeloffe'  (1983)  explanation  of  a 
timing  error  is  plausible.  However,  wide-angle  re¬ 
flection  measurements  have  been  made  with  radar 
sounding  systems  other  than  the  ones  they  state  to 
have  the  timing  error  and  still  a  discrepancy  exists  in 
the  D  and  e'  determinations.  We  suggest  that 
curved  ray  path  effects  associated  with  Snell's  law 
and  the  standard  wide-angle  reflection  analysis  us¬ 
ing  a  time2-distance2  plot  are  a  cause  of  the  discrep¬ 
ancy.  Indeed,  the  findings  of  Hughes  (1985)  and 
Morey  and  Kovacs  (1985)  show  that  a  standard 
wide-angle  reflection  measurement  analysis  gives 
lower  e'  and  higher  D  determinations  than  are  ap¬ 
propriate. 
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APPENDIX  A:  SUPPLEMENT  EQUATIONS  FOR  DETERMINING 
THE  DIELECTRIC  CONSTANT  OF  AN  AIR-ICE  MIXTURE 


Additional  formulations  are  presented  in  this  section  for  determining  the  e'  vs.  p  for  fim 
and  ice.  All  equations  have  been  offered  by  the  indicated  authors  as  a  means  of  determining  the 
e' of  an  air-ice  mixture.  Below  each  equation  is  a  table  giving  example  air-ice  mixture  specific 
gravities  and  the  calculated  e'r  value  as  determined  by  eq  27  and  by  the  new  expressions  pre¬ 
sented. 

Tiuriet  al.  (1984)  made  measurements  in  the  0.85-12.6  GHz  range  on  snow.  Their  empirical 
equation  for  e'  is 


e' =  1  +  1.7  p  +  0.7  p2 
which  they  state  for  practical  purposes  can  be  simplified  to 

e;=l  +  2p. 


(Al) 


(A2) 


p 

Ef  (eq  27) 

er  (eqAl) 

e;(eqA2) 

0.0 

1.0 

1.0 

1.0 

0.2 

1.37 

1.37 

1.40 

0.4 

1.79 

1.79 

1.80 

0.6 

2  27 

2  27 

2.20 

0.8 

2.81 

2.81 

2.60 

0.917 

3.15 

3.15 

2.83 

Ambach  and  Denoth  (1980)  give  the  following  empirical  equation  for  their  dry  snow 
measurements  at  20  MHz: 


£'  =  1  +  2.2  p. 


(A3) 


P 

er  (eq  27) 

e;(eqA3) 

0.0 

1.0 

1.0 

0.2 

1.37 

1.44 

0.4 

1.79 

1.88 

0.6 

2.27 

232 

0.8 

2.81 

2.76 

0.917 

3.15 

3.02 

Hallikainen  et  al.  (1982)  made  measurements  at  4-18  GHz  on  dry  snow  and  found  empiri¬ 
cally  that 


e;=  1.91  p. 


(A4) 
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p 

e;<eq27) 

e;(eqA4) 

0.0 

IB 

1.1 

02 

1 37 

138 

0.4 

1.79 

1.76 

0.6 

227 

2.15 

0.8 

2.81 

253 

0.917 

3.15 

2.75 

Bums  et  al.  (1985)  made  measurements  at  100  MHz  on  diy  snow  and  found  empirically  that 


e'r=  1.1  +  2.2  p. 


(A5) 


P 

e;(eq27) 

e;(eqA5) 

0.0 

1.0 

1.1 

02 

137 

154 

0.4 

1.79 

1.98 

0.6 

227 

2.42 

0.8 

2.81 

256 

0.917 

3.15 

3.12 

Fujita  et  al.  (1992)  made  measurements  at  9.7  GHz  on  polycrystalline  ice  and  found  from 
their  results  that 


ej.  =  3.08p  +  0.41. 


(A6) 


P 

e;(eq27) 

e;(eqA6) 

0.0 

1.0 

0.41 

02 

137 

1.03 

0.4 

1.79 

1.64 

0.6 

227 

226 

0.8 

2.81 

2.87 

0.917 

3.15 

323 

Pearce  and  Walker  (1967)  proposed  the  following  expression  based  on  theoretical  con¬ 
siderations: 


e'  =  3.16  p  +  0.41 


(A  7) 


where  0.535  2  p  £  0.920. 
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p 

e;(eq27) 

e;(eqA7) 

2 27 

231 

0.7 

253 

2.62 

0.8 

251 

2.94 

3.15 

331 

Sihvola  et  aL  (1985)  proposed  a  mixing  formula  by  Taylor  (1965)  for  an  air-ice  mixture  with 
randomly  oriented  disk-shaped  inclusions  of  the  form 


♦  -V  (» -  *«) 

3eri 


which  can  be  simplified  to 


e'  =  (1.007  +  0.838 p)2  (A9) 

with  a  cross-correlation  of  r2 = 1 .000.  Sihvola  et  aL  showed  tha  t  eq  A8  fitted  their  data,  obtained 
at  0.8-13  GHz,  better  than  another  one  by  Taylor  for  spherical  inclusions. 


P 

£'r(eq27) 

e;(eqA9) 

0.0 

1.00 

1.01 

02 

137 

138 

0.4 

1.79 

1.80 

0.6 

227 

228 

03 

231 

231 

0.917 

3.15 

3.15 

In  an  important  paper  on  dielectric  mixing  formulas,  Sihvola  (1989)  presents  a  self-consis¬ 
tent  general  mixing  formula  for  determining  the  complex  permittivity  of  a  homogeneous  ma¬ 
terial  containing  spherical  inclusions.  As  applied  to  determining  the  dielectric  constant  of  dry 
fim  and  ice,  his  expression  is 


Er  = 


ib2  +  4oc  -  b 


2a 


(A10) 


where 


b  =  e'r  j  +  2ea  -  2oea  -/*  (e'ri  -  ea)  (1  +  o) 
c  =  ea[e'i  +  (2  -  o)ea  +/s  (2  -  a)(e'ri  -  ea)] 
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and  parameter  o  is  a  positive  integer  and  parameter /;  is  the  ice  density  fraction. 

Sihvola  showed  the  mterrelationship  between  his  mixing  formula  and  that  of  other  authors. 
His  results  indicated  that  when  o  =  3,  eq  A10  gave  e'  values  for  fim  and  ice  in  close  agreement 
with  eq  Al.  Equation  A10,  with  0  =  3,  can  be  simplified  to 

e;  =  (0.988  +  0.859p)2  (All) 

with  a  cross-correlation  of  r2  =  1.000. 


p 

e'r(eq27) 

er(eqAll) 

0.0 

1.00 

0.98 

02 

137 

135 

0.4 

1.79 

1.77 

0.6 

227 

226 

0.8 

2.81 

2.81 

0.917 

3.15 

3.15 

Sihvola  also  showed  that  an  expression  by  Sen  et  al.  (1981)  for  a  medium  containing  spheri¬ 
cal  inclusions: 


(A12) 


when  applied  to  fim  and  ice  gave  e'  values  in  good  agreement  with  equation  Al.  When  eq  A12 
is  simplified  to  a  refraction  type  expression  we  obtain 


e;  =  (0.995  +  0.848 p)2 
with  a  cross-correlation  of  r2  =  1.000. 


(A13) 


P 

e;(eq27) 

er(eqA13) 

0.0 

1.0 

0.99 

02 

137 

136 

0.4 

1.79 

1.78 

0.6 

227 

226 

0.8 

2.81 

2.80 

0.917 

3.15 

3.14 

The  interested  reader  is  directed  to  Sihvola's  paper  for  additional  mixing  formulas  that  may 
be  applied  to  estimating  the  dielectric  constant  of  fim  and  ice. 

Of  foe  above  expressions,  it  appears  that  eq  Al,  A9,  All  and  A13  correlate  best  with  eq  27. 
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APPENDIX  B:  THE  DIELECTRIC  CONSTANT  OF  A  SOIL-WATER  MIXTURE. 


The  dielectric  constant  of  dry  soil  e'^  typically  lies  between  2  and  4,  while  the  relative 
dielectric  constant  of  water  e'  w  is  81  at  20°C.  This  value  will  vary  with  temperature  from  about 
88  at  0°C  to  lower  values  as  the  temperature  climbs.  The  value  for  e' w  is  also  frequency  depen¬ 
dent.  Since  e'wis  such  a  large  value  it  has  a  major  influence  on  the  relative  dielectric  constant  of 
the  soil-moisture  mixture  e'm.  Here  we  show  the  variation  of  e' m  vs.  soil  moisture  content  mv 
at  a  fixed  frequency  of  1.074  GHz  and  at  room  temperature.  The  soil  moisture  content  is  the 
ratio  of  the  weight  of  water  per  unit  volume  of  soil. 

Presented  in  Figure  B1  are  the  results  of  403  laboratory  e'm  vs.  mv  tests  by  Lundien  (1971). 
His  measurements  were  made  on  sand,  silt,  silt-loam  and  clay  soils  using  an  L-band  interfer¬ 
ometer.  Also  shown  is  a  refraction  equation  which  well  describes  e' w  vs.  mv  for  the  com¬ 
mingled  soils.  The  regression  curve  in  Figure  B1  fits  the  data  with  an  r2  =  0.970  and  a  standard 
error  of  ±  2.02.  When  the  regression  curve  is  forced  through  an  e'  w  for  water  of  81,  the  equation 
for  the  curve  becomes  as  shown  in  Figure  B2.  The  data  lying  below  the  regression  curve  tends 
to  be  those  for  day  type  soil  while  the  data  above  foe  curve  is  that  for  foe  coarser  grained  soils. 
Representative  regression  curves  passed  through  foe  sand  and  clay  material  and  forced 
through  an  e'  w  of  81  at  mv  =  1  are  shown  in  Figure  B3. 

The  test  results  shown  in  Figure  B1  will  vary  with  temperature  as  previously  mentioned. 
However,  Lundien  (1971)  found  that  over  a  temperature  range  of  0°  to  65°C  foe  change  in  e'm 
appears  to  be  uniform,  varying  at  about  1*2%  per  °C. 


Figure  BI.  Dielectric  constant  of  soil  vs.  moisture  content. 


21 


Figure  B2.  Dielectric  constant  of  soil  vs.  moisture  content  where 
the  regression  curve  was  forced  through  a  dielectric  constant  for 
water  of  81. 


Figure  B3.  Representative  curves  for  the  dielectric  constant  of 
sand  and  clay  soil  vs.  moisture  where  the  respective  curves  were 
forced  through  a  dielectric  constant  for  water  of  81. 
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e'  values  with  Robin's  equation,  laboratory  measurements  on  fim  and  ice  and  other  expressions  given  in  fire  litera¬ 
ture  for  determining  e'  vs.  the  specific  gravity  of  dry  fim  and  ice.  Our  findings  indicate  that  the  form  of  Robin's  equa¬ 
tion  is  valid.  However,  our  analysis  also  indicates  the  expression  could  be  slightly  improved  to  read  e'  =  (1  +  0.845 
p)2.  Reasons  are  suggested  as  to  why  previous  wide-angle  radar  sounding  studies  did  not  reproduce  Robin's  find¬ 
ings. 
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Throughout  the  report,  £'  should  read  £e  and  be  defined  as  the  effective  (bulk)  dielectric  constant. 
Inside  cover 

Abstract,  line  5,  "1968"  should  read  "1969" 
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